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In the Dirac theory for the motion of free relativistic electrons, highly oscillatory components 
appear in the time evolution of physical observables such as position, velocity, and spin angular 
momentum. This effect is known as zitterbewegung. We present a theoretical analysis of rather 
different Hamiltonians with gapped and/or spin-split energy spectrum (including the Rashba, Lut- 
tinger, and Kane Hamiltonians) that exhibit analogs of zitterbewegung as a common feature. We 
find that the amplitude of oscillations of the Heisenberg velocity operator v(t) generally equals the 
uncertainty for a simultaneous measurement of two linearly independent components of v. It is 
also shown that many features of zitterbewegung are shared by the simple and well-known Landau 
Hamiltonian describing the dynamics of two-dimensional (2D) electron systems in the presence of a 
magnetic field perpendicular to the plane. Finally, we also discuss the oscillatory dynamics of 2D 
electrons arising from the interplay of Rashba spin splitting and a perpendicular magnetic field. 

PACS numbers: 73.21.-b, 71.70.Ej, 03.65.Pm 
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I. INTRODUCTION AND OVERVIEW 

The Dirac equation^^ was derived to obtain a rela- 
tivistic generalization of Schrodinger's approach to quan- 
tum physics that describes the dynamics of single- 
electron quantum states. While it served as an impor- 
tant stepping stone towards a more complete description 
of quantum-electrodynamic effects, Dirac theory has oc- 
casionally been regarded with some suspicion. In partic- 
ular, the effect of zitterbewegung 4 showed that solutions 
of the Dirac equation exhibit peculiarities that are in- 
consistent with classical intuition in a more fundamental 
way than nonrelativistic quantum physics. The zitter- 
bewegung is an oscillatory dynamics of observables in- 
duced by the Dirac equation, with a frequency of the 
order of 2mc 2 /%, where m is the electron mass, c is the 
speed of light and Ti is the Planck constant. The ampli- 
tude of oscillations in a particle's position is of the order 
of the Compton wave length. Subsequently, zitterbewe- 
gung attracted some interest as a possible way to under- 
stand the intrinsic magnetic moment of the electron^ 
Later, on the level of fundamental physics, the advent 
of quantum field theory obviated the need to discuss rel- 
ativistic quantum theory in terms of a first-quantized, 
Schrodinger-type theory. Present interest in the Dirac 
equation ranges from hadronic physics^ over lattice gauge 
theory^ to recent efforts^ to incorporate relativistic ef- 
fects into quantum-chemistry calculations. 



A Dirac-like dynamics causing analogs of zitterbewe- 
gung was also predicted for electrons moving in crys- 
talline solids ) 10 : 11 in particular for narrow-gap semi- 
conductors^ carbon nanotubes,— graphene sheets^ 
tunnel-coupled electron-hole bilayersi^ and superconduc- 
tors^ All these systems are characterized by having the 
relevant electron excitations grouped into two bands sep- 
arated by a nonzero energy gap so that their energy spec- 
trum is similar to the spectrum of the Dirac Hamiltonian. 
A recent studj*il of two-dimensional (2D) electron sys- 
tems in inversion-asymmetric semiconductor heterostruc- 
tures showed the presence of an oscillatory motion anal- 
ogous to zitterbewegung arising from spin splitting of the 
energy levels. The spin splitting corresponds to an en- 
ergy gap that vanishes for momentum p — > 0. A similar 
situation occurs for electronic excitations in the bulk of 
an ideal graphene sheet ^ 

These findings indicate the need to understand zitter- 
bewegung-like effects on a more general level. In Ref. [Til , 
the authors presented a general formula for the Heisen- 
berg position operator r(t) in systems that can be de- 
scribed by effective 2x2 Hamiltonians ^ In the present 
work, we have investigated the oscillatory dynamics of 
Heisenberg observables such as position r(t), velocity 
v(i) = dv/dt, orbital angular momentum L(i), and spin 
S(t) in a variety of qualitatively different models that 
describe the motion of free (quasi-)particles. Besides 
the Dirac Hamiltonian, we have studied three Hamiltoni- 
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ans frequently used in semiconductor physics to describe 
the dynamics of (quasi-free) Bloch electrons in the vicin- 
ity of the fundamental gap, the Rashba, 20 Luttinger, 21 
and Kane 22 Hamiltonians. A number of striking features 
emerge quite generally in all these models, thus illustrat- 
ing remarkable similarities between time evolutions gen- 
erated by rather different Hamiltonians. We suggest that 
these common features can be used to extend the concept 
of zitterbewegung to a broader class of quantum Hamilto- 
nians for free (quasi-)particles. Our analysis shows that 
this generalized notion of zitterbewegung is manifested, 
in addition to the oscillatory unitary time evolution of 
observables, also by uncertainty relations characterizing 
the measurement of such observables. These two aspects 
turn out to be closely related. In particular, they can 
be described, for each of the models considered here, 
by the same set of parameters. Also, we identify the 
typical scales (lengths, velocities, and frequencies) that 
characterize zitterbewegung-like oscillatory motion. We 
emphasize that this extended notion of zitterbewegung is 
entirely based on quantum mechanical concepts. In an 
alternative, semiclassical approach one would identify a 
zitterbewegung relative to a suitable classical dynamics 
as, e.g., in Ref. [23]. In some cases the conclusions will 
be different from those obtained within the present ap- 
proach. The most general aspects of our study can be 
summarized as follows: 

(i) An oscillatory motion occurs in the time evolution 
of free (quasi-) particles when the energy spectrum of the 
corresponding Hamiltonian H is characterized by one or 
several energy gaps. Besides the Dirac model, an im- 
portant example are Bloch electrons in solids, 10 ' 11 whose 
quantum dynamics are described by effective free-particle 
Hamiltonians that incorporate the effect of the periodic 
lattice potential. 

(ii) In the case of two-band models (e.g., the Dirac, 
Rashba, and Luttinger models), zitterbewegung-\\ks ef- 
fects are generally characterized by an amplitude op- 
erator F and a frequency operator d>(p). These two 
quantities enter the expression for the velocity operator 
in the Heisenberg picture, which can be decomposed as 
v(t) = v(i) + v(i), where the mean part is 



r(t) = r(t) + r(t), where 



. . dH „ 
= ~dp ~ 

and the oscillating part is 

v(t) — F e~ il2 '( p " ) * = e^P) * F 



(la) 



(lb) 



Here, hu)(p) is related to the energy difference between 
states having the same momentum p, but belonging to 
different subspaces (i.e., energy bands) of the Hamilto- 
nian. The operator F, which anticommutes with c£>(p), 
determines the magnitude of oscillations in the velocity 
components but also enters the expression for the mean 
part. We can integrate Eq. |T]) to get the Heisenberg po- 
sition operator that can be decomposed in the same way, 



r(i) = r + vi + F 



f(t) = F 



,— iui(p)t 

iu>(p) 



Similarly, we get the time derivative of v(t), 
v(i) = ia>(p)F e-^P)* . 



(2a) 
(2b) 

(3) 



The operators F and u>(p) govern also the oscillations 
in the Heisenberg time evolution of the orbital angular 
momentum operator L(t), and the spin operators S(t). In 
systems with more than two bands (Kane and Landau- 
Rashba models), more than one characteristic frequency 
and amplitude operator can appear. 

(iii) For each model describing an oscillatory multiband 
dynamics of free particles, the components of the velocity 
operator v(i) do not commute. This can be written as 
an uncertainty relation that takes the form (apart from 
a prefactor of order one) 



Avj Av k > w 2 (j ^ fc), 



(4) 



where v is the amplitude of the oscillatory motion, see 
Eq. (|lbp. The uncertainty relations ^ are an integral 
part of our analysis^ 

(iv) The velocity operator v(t) does not commute with 
the Hamiltonian. Although we discuss the motion of free 
(quasi-)particles, the components of v(t) are not con- 
stants of the motion, see Eq. ([3]). On the other hand, 
momentum p is always a constant of the motion. This 
implies that none of the models discussed here provides 
a simple relation between momentum p and velocity v. 

(v) The counterintuitive properties of r(t) and v(t) 
arise because r(t) mixes different subspaces Ttj that are 
associated with the different bands in the energy spec- 
trum of H. Thus we can interpret zitterbewegung-likc 
phenomena as an interference effect. In the case of two- 
band models, one can replace r by the part r that leaves 
the subspaces 7i± associated with the '+' and '— ' bands 
separately invariant, 



(5) 



where P± are projection operators onto these subspaces. 
The result coincides with the mean part r(t) of r(t) in- 
troduced in Eq. (|2a,|) . i.e., zitterbewegung-likc effects arc 
removed by the projection ([5]). This result can be un- 
derstood from a different perspective by analyzing the 
amplitude operator F. We get 



FP+ = P_F and FP_ = P+F , 



(6) 



i.e., F maps states associated with the '+' band onto 
states associated with the '— ' band and vice versa. An 
alternative definition of r(t) is obtained by applying the 
inverse unitary transformation to r(t) that makes H di- 
agonal. The same techniques can also be applied to v(i) 
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to obtain v(t) given in Eq. (fla[) . The components of v(t) 
commute; hence they can be measured simultaneously 
[unlike Eq. J3J]. They also commute with the Hamilto- 
nian so that they are constants of the motion. 

(vi) In every case considered, zitterbewegung-like phe- 
nomena are manifested also by oscillations of the orbital 
angular momentum L(t) and spin S(t). At the same time, 
the total angular momentum J does not oscillate as a 
function of time. As expected for a model of a free parti- 
cle, J is a constant of the motion, i.e., it commutes with 
the Hamiltonian. From a different perspective, this im- 
plies that the oscillations of L(t) = r(t) X p and S(t) 
must cancel each other, which is possible only if the os- 
cillations of r(t) and S(t) have a common origin. For 
the Rashba Hamiltonian, the oscillatory motion of S z (t) 
corresponds to the well-known and experimentally ob- 
served^ spin precession in the effective magnetic field of 
the Rashba term. 

The following Sections HIHVl are devoted to a de- 
tailed discussion of zitterbewegung effects arising in sys- 
tems whose time evolution is governed by the Dirac, 3 
Rashbaj2£ Luttinger^i and Kane 2 ^ Hamiltonians. Re- 
markable formal similarities between the oscillatory be- 
havior of observables in these models are established, as 
outlined above. Next we show in Sec. IVII that the fa- 
miliar Landau model of 2D electrons subject to a per- 
pendicular magnetic field 2 ^ exhibits essentially all the 
features attributed to the extended notion of zitterbewe- 
gung in previous sections. We finish our case studies in 
Sec. lVIll bv investigating the quantum oscillatory dynam- 
ics of 2D electrons arising from the interplay of Rashba 
spin splitting and a perpendicular magnetic field^ Con- 
clusions and a summary of open questions are presented 
in Sec. IVIIII For easy reference, we provide a number of 
relevant basic formulae in the Appendix. 



Note that H D 
velocity operator reads 



h 2 D . In the Schrodinger picture, the 



dH D 
dp 



(10) 



so that the components of v have the two discrete eigen- 
values ±c. In the Heisenberg picture, we get Eq. with 



dHr 



2 

C p 



- TJ— , £(P) 



2H T 



dp E D 



(11) 



The operator F mediates a coupling between states with 
positive and negative energies; see below. The oscilla- 
tory part v(t) of v(t), given in Eq. (fTb|) . describes the 
zitterbewegung. The frequency of the zitterbewegung is 
(at least) of the order of ujd = 2mc 2 /h. Integrating v(i) 
yields the position operator r(t) in the Heisenberg pic- 
ture, see Eq. ([2]), which contains again the quickly oscil- 
lating term e~ lw ( p )*. The oscillatory time dependence is 
similar to the motion of a nonrelativistic particle in the 
presence of a magnetic field [see Eq. (|4"5]l and discussion 
in Sec. IVIj . An illuminating discussion of zitterbewegung 
based on a numerical calculation of the time evolution of 
wave packets can be found in Ref. |29j. 

It turns out 3 - that orbital angular momentum L = r x p 
and spin S, which is defined as 



ifi h 

S = a x a = — 

4 2 



a 
a 



(12) 



show the phenomenon of zitterbewegung, too. For the 
orbital angular momentum, we have 



L(t) = r(i) xp = L + Fxp 



1 _ e -iii(p)t 



(13a) 



The time evolution of spin in the Heisenberg picture reads 



II. DIRAC HAMILTONIAN 

Our discussion of zitterbewegung for the Dirac Hamil- 
tonian follows, for the most part, Ref. 3. We include this 
section with an overview of the effect's salient features to 
provide a reference frame and notation for our following 
discussion of solid-state analogies. 

The Dirac Hamiltonian Ho describes a free relativistic 
electron or positron. It can be written in the form 



Ho = col ■ p + P mc 



where 



a = 



a 
a 



(3 = 



1 2 X2 
-1 2 X2 



(7) 



(8) 



and cr is the vector of Pauli spin matrices. (Here we 
assume magnetic field B = 0. See Ref. [28| for the gener- 
aliziation to finite B.) The energy eigenvalues of Ho are 
-^±(p) — ±Ad, where 



A 



D 



c 2 p 2 . 



(9) 



S(t) = S — P x p 



1 



>( P )t 



(13b) 



Thus it follows from Eqs. (|13ap and (|13b[) that the total 
angular momentum J = L + S does not oscillate as a 
function of time, 



J(i) = J = L + S 



(13c) 



which reflects the fact that [J, Ho] = 0. 

We can estimate the magnitude of zitterbewegung by 
evaluating the square of v(t) (Ref. |j)- This yields 



2 c 2 (2A 2 3 +m 2 c 4 ) 



*(*) = 



A 2 

D 



(14) 



i.e., v 2 varies between 3c 2 in the nonrelativistic limit and 
2c 2 in the relativistic limit. (Note that, although v 2 > c 2 , 
no measurable velocity exceeds c.) On the other hand, 
the components of the velocity operator v do not com- 
mute. Equations (JTUJ) and (fl"2")l imply that 



r l 4iC o 

[Vj,Vk\ = -r— £jkl £>l 



(15a) 
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Diagonalizing this equation yields the uncertainty rela- 
tion for j =/= k 



Av 3 Av k > c 2 = (i LO D X D ) 



,2 



(15b) 



where = h/(mc) is the Compton wave length. Thus 
both the magnitude and the uncertainty of the zitterbe- 
wegung are given by c 2 (Ref. 24). We can also estimate 
the spatial amplitude of the zitterbewegung using the de- 
composition r(t) = f (t) + f(t), see Eq. ([2]). We get for 
the oscillating part 



r 2 



(*) 



h 2 c 2 



(2A 2 D 



m 2 c ) 



4A% 



(16) 



i.e., in the nonrelativistic limit, the amplitude of zitter- 
bewegung is approximately Xd, and it is given by the de 
Broglie wave length Xb — fo/p in the relativistic limit. 

It is well-known^ that zitterbewegung is caused by a 
coupling between the states with positive energies ( "par- 
ticles", subspace TL+) and negative energies ("antiparti- 
cles", subspace TC-). Thus one can eliminate the oscilla- 
tions of r(t) by projecting r on Ti.± as in Eq. ([5]) and the 



result coincides with Eq. 
do not commute: 



a|. The components r*j of r 



he 2 



[rj,r k \ 



A D 



where 



S F x p 



c2)(p) 



(17a) 



(17b) 



is the spin operator analogous to Eq. (fTa[) that does not 
mix the subspaces of positive and negative energy states. 
Diagonalizing Eq. (|17ap yields the uncertainty relation 
for j =/= k 



Afj Af fc > 



Ti 2 c 2 yjm 2 c± + \e ]kl \c 2 p 2 



A3 
1V D 



(17c) 



Thus we obtain in the nonrelativistic limit ( li v <C c" ) 

Af Af k >\X 2 D . (17d) 

In the opposite (ultrarelativistic) limit ( "v ^ c" ) we have 
AfjAf fc >iA| . (17e) 

The time derivative v of the mean position operator r 
is the velocity operator one would expect based on the 
correspondence principle and classical relativistic kine- 
matics. Its components Vj, v k commute, [vj , Vfc] = 0. 
Furthermore, v commutes with Hp, i.e., it is a constant 
of the motion. 

Equation |5j) is motivated by the requirement that 
it leaves the subspaces TL± separately invariant. How- 
ever, this requirement is not sufficient for a unique defini- 
tion of a relativistic position operator. The Dirac equa- 
tion becomes diagonal in the Foldy-Wouthuysen (FW) 
representation^ If we require the mean position opera- 
tor to be diagonal in this representation, it can be ob- 
tained in the standard representation via an inverse FW 
transform^ 



r NW (*) = r + vt 



hp 

2iA 



D 



COL — 



c 3 (a. ■ p) p 
A D (A D + m& 

I 



c 2 S x p 



A D (A D +mc 2 ) 



(18) 



This operator is often called the Newton- Wigner position 
operator^ In contrast with the position operator r in 
Eq. ([5]), the components of r NW do commute. We see 
that f NW (t) shares with r(t) from Eq. © that the time 
derivative v(t) is given by Eq. ([la)), i.e., it does not show 
any zitterbewegung because it is a constant of the motion. 



General requirements for any position observable de- 
scribing the localization of a particle or wave packet are 
discussed in Ref. [||. In this context, the operator r in 
Eq. (|5|) appears inappropriate because its components do 
not commute. The optimal choice for a position observ- 
able is the operator r NW . However, general arguments 
prohibit the possibility of strict spatial localization for 
a one-particle state (see, e.g., Refs. 3.32.33). This im- 
poses restrictions on the utility of any position operator 
in relativistic systems. 



III. RASHBA (AND PAULI) HAMILTONIAN 

An intriguing example of zitterbewegung-\ike dynam- 
ics exhibited by a non-Dirac-like Hamiltonian has been 
founcUi in the Rashba model^S, This model describes 2D 
electrons in semiconductor heterostructures with spin- 
orbit coupling present, using the effective Hamiltonian 
(we assume here B = 0, see Sec. lVIII for the case B ^ 0) 



H = 



P„ 



2m 



Here 



H R = a (a- x p) 



e z = a 



H, 







Py - IP* 



Py + iPx 




(19a) 



(19b) 



is the Rashba term with Rashba coefficient a (Ref. \2w , 
and e z denotes the unit vector in the direction perpen- 
dicular to the 2D plane. (Note that H\ — a 2 p 2 , similar 
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to the Dirac Hamiltonian.) The Hamiltonian (|19|) is also 
equivalent to the Pauli Hamiltonian^ for a 2D system. 
The energy eigenvalues of H are 



E±(p) 



P 

2m 



± ap 



(20) 



The time-dependent position operator r(i) in the 
Rashba model was discussed previously in Ref.ll7l. Evalu- 
ated in close analogy to the Dirac case, it is again possible 
to decompose r(t) into a mean part f (i) and an oscillat- 
ing part f (t). The result is of the form shown in Eq. |2]), 
where F and cD(p) are now 



dp 



2 

crp 



a 2 e z 



x P 



Hr 



L){p) 



2H R 



(21) 

Explicit evaluation shows that r(t) oscillates with the 
frequency ojr — 2ap/h, which is equal to the precession 
frequency of a spin moving in the effective magnetic field 
of the Rashba term [see Eq. (|24b[) below] . The oscillation 
becomes arbitrarily slow for p — > 0. We find for the 
oscillating part of r(i) 



f 2 (t) = X 2 B /4 



(22) 



i.e., the magnitude of the oscillations is of the order of 
the de Broglie wave length A^ and independent of the 
Rashba coefficient a. Note that As diverges in the limit 
p^O. 

We obtain the mean part f(t) by projecting on the 
subspaces of H associated with the spin-split bands, as 
in Eq. ([5]). We find the same r(t) by applying an inverse 
FW transformation, similar to Eq. (TT5)) . For the Rashba 
model, the last term in Eq. (|2"aj) corresponds to a spatial 
separation of up and down spin contributions in a wave 
packet by ~ As (independent of the Rashba coefficient 
a), which was noticed in previous numerical work^ The 
general validity of Eq. (f2"aj) for two-band models implies 
the existence of similar displacements for the Dirac and 
Luttinger cases. See also Ref. [l8|. The components x 
and y of the mean position operator r commute, similar 
to r NW in Eq. (| 18[) . 

The velocity operator and its derivative are given by 
Eqs. H]) and ©, respectively, using expressions (|2"Tj) . The 
oscillatory part of v satisfies 



v\t) 



a 



i.e., the magnitude of the oscillatory motion v(i) is given 
by the Rashba coefficient a. On the other hand, the 
components of v do not commute, and we have 

[v x {t),v v (t)}=2ic?a z e- i *^ t , (23b) 

which implies 

Av x Av v >a 2 = {\u R X B f , (23c) 

analogous to Eqs. (JTSJ) - In Eq. (|23cp we replaced the 
matrix-valued RHS of Eq. (|23b[) by the eigenvalues of 



this matrix. Thus similar to the Dirac case, both the 
magnitude of the oscillations in v(£) and the minimum 
uncertainty are given by the same parameter. The com- 
ponents of the mean part of the velocity operator com- 
mute, [v x (t),v y (t)] = 0. They also commute with Hr, 
i.e., they are constants of the motion. 

The time dependence of orbital angular momentum L z , 
spin component S z , and total angular momentum J z = 
L z + S z can be straightforwardly discussed. We get 



+ *Zl(l- e-^(p)*) , (24a) 



2 



-iuj(p)t 



s z 



(24b) 
(24c) 



The formal structure of these equations is analogous to 
the Dirac-case counterparts shown in Eqs. (fTS"]) . Equa- 
tion (|24b[) represents the well-known spin precession in 
the effective magnetic field of the Rashba term, which has 
been observed experimentally.— The total angular mo- 
mentum component perpendicular to the plane does not 
depend on time, as expected from [J z , H] = 0. Obviously 
Eqs. (|24|) require that the spin precession is caused by the 
effective in-plane magnetic field of a spin-orbit coupling 
term such as the Rashba term. We see here clearly the 
difference between spin precession caused by spin-orbit 
coupling and spin precession caused by the Zeeman term 
in the presence of an external in-plane magnetic field. In 
the latter case J z is not a constant of motion. 



IV. LUTTINGER HAMILTONIAN 

The uppermost valence band Tg of common semicon- 
ductors like Ge and GaAs is well-characterized by the 
Luttinger Hamiltonian^ 



H 



7i P 
2m 



Hr 



(25a) 



We assume B = and use the spherical approximatio: 



,35 



Hr = 



1_ 

m 



(p-S) 2 -|p 2 l 4 x4 



(25b) 



where 71 and 7 are the dimensionless Luttinger param- 
eters, and S is the vector of 4 x 4 spin matrices for a 
(23a) system with spin s = 3/2. (Note H\ = ; y 2 p 4 /m 2 , sim- 



ilar to the Dirac Hamiltonian.) The twofold-degenerate 
energy eigenvalues of H are 



£±(p) = -f-(7i±2 7 ) 
Im 



(26) 



The upper sign corresponds to the so-called light-hole 
(LH) states with spin-z component M = ±1/2, and the 
lower sign corresponds to the heavy-hole (HH) states with 
M = ±3/2. The momentum-dependent energy gap be- 
tween HH and LH states is 



Tiuji 



2~fp 2 jra 



(27) 
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The position operator is of the form shown in Eq. (|2j) 
with F and lu(p) given by 



dp 



2pH L 

p2 



Lb(p) 



2H L 



(28) 



Thus r(i) oscillates with the frequency ojl, which has 
been noticed in previous numerical worki 3 ^ Similar to 
the Rashba Hamiltonian, these oscillations become arbi- 
trarily slow for p — > 0. The squared amplitude of the 
oscillations of r(t) is r 2 (t) = (3/2) A^, independent of the 
Luttinger parameter 7. It diverges for p — > 0. 

We obtain the mean position operator r, defined in 
Eq. ([5]), using projection operators that project onto HH 
and LH states^ The result coincides with Eq. (|2"a|). The 
components of f do not commute, 



[rj,r k \ 



hdH L /d Pj hdH L /d Pk 



2Hr 



2Hr 



implying the uncertainty relation 



Afi Affc > 



3_£jfcj h pi 
4p3 



(29) 



(30) 



This uncertainty is of the order of (or less than) the 
de Broglie wave length. The uncertainty is the largest 
for those components fj that are perpendicular to p. 



Using Eqs. (fzgj) . the velocity operator can be written 
in the form shown in Eq. |T]). For its oscillating part, we 
find v 2 (t) = 6j 2 (p/m) 2 . The components of v do not 
commute, 



[Vj,Vkl 



OHr OHr 



dpj ' dp k 



(31a) 



which corresponds to the uncertainty relation 



Av x Av y > 



1 2 P 



3(4p| 



4p2 



and cyclic permutations thereof, i.e., the uncertainty is 
approximately limited by 37 2 (p/m) 2 = j(ujl^b) 2 - Thus 
again, the magnitude of the oscillations of v(i) and the 
minimum uncertainty are characterized by the same com- 
bination of parameters. The velocity v(t) is not a con- 
served quantity but satisfies Eq. (J3j) . However, the mean 
velocity operator is again given by Eq. {Taj. Its compo- 
nents commute and are constants of the motion. 

The time dependence of orbital angular momentum L, 
spin S, and total angular momentum J = L + S turns 
out to be given by Eqs. (fT3")) . Note that the time depen- 
dence of S(t) in the Luttinger model corresponds to a 
spin precession in the absence of any external or effec- 
tive magnetic fields Again, the total angular momen- 
tum does not depend on time, which reflects the fact 
that [J,H] = 0. 

We remark that a similar analysis as presented in this 
Section also applies to models that neglect the spin de- 
gree of freedom. An example is the 3x3 Shockley Hamil- 
tonian that describes spinless holes in the uppermost va- 
lence band Tg of semiconductors like Si. 39 ' 40 Indeed, this 



is consistent with the fact that zitterbewegung for the 
Dirac case can be studied already in a model with only 
one spatial dimension, where the Dirac Hamiltonian Hrj 
becomes a 2 x 2 matrix that reflects the occurence of 
both signs of the energy in the spectrum of Ho', but 
this Hamiltonian does not describe the spin degree of 
freedom.- 29 A spin with spin-orbit coupling is not a nec- 
essary condition for the oscillatory behavior of r(t) and 
v(i) to occur. The most basic ingredient required for zit- 
terbewegung-\ike effects are several bands separated by a 
(usually momentum-dependent) gap. Often the splitting 
of these bands can be described by an effective spin-orbit 
coupling.— i 41 - 



V. KANE HAMILTONIAN 

The Kane Hamiltonian^ 2 - is an effective Hamiltonian 
that captures the important physics of electrons and 
holes in narrow-gap semiconductors like InSb. We re- 
strict ourselves to the 6x6 Kane model which includes 
the lowest conduction band Tg and the uppermost va- 
lence band Tg (3D, B = 0), neglecting the split-off va- 
lence band r*, because that model permits a fully ana- 
lytical solution. Then we have 



H 



K 



(E g /2) l s 



P 



V3PT-P 

-(E g /2) 1 4 X4 



(32) 



Here E g is the fundamental energy gap, and V denotes 
Kane's momentum matrix element. The vector T of 2 x 4 
matrices is defined in Ref. l42l The energy eigenvalues of 
Hk are (each twofold degenerate) 



E±(p) = ±A K , E {p) = -E g /2, 



(31b) where 



A 



K 



^(E g /2f 



\V 2 p 2 



(33a) 



(33b) 



i.e., the Kane Hamiltonian combines the gapped spec- 
trum of the Dirac Hamiltonian with the gapless spectrum 
of the Luttinger Hamiltonian. (Indeed, the Luttinger 
Hamiltonian corresponds to the limiting case E g — > 00 of 
the Kane model.) The energy spectrum E±^o(p) is shown 
in Fig. m 

Results similar to those discussed below can also be 
derived perturbatively for the full 8x8 Kane Hamilto- 
nian that includes the split-off valence band Ti; . We also 
remark that a simplified 4x4 Kane Hamiltonian, which 
includes only the conduction band Tg and the valence 
band Ty, is strictly equivalent to the Dirac Hamiltonian. 
Recently, zitterbewegung was studied for a simplified ver- 
sion of the Kane model where the HH band [with dis- 
persion Eq(p) = —Eg/2) and the split-off band were 
neglected^ In this limit, the Kane Hamiltonian becomes 
similar to the Dirac Hamiltonian. Our analysis below 
shows that qualitatively new aspects arise when the HH 
band is taken into account. 




(e) 




FIG. 1: Energy spectrum £±,o(p) of the 6x6 Kane model. 
Here, E g denotes the fundamental gap. Each band E(p) is 
twofold degenerate. 



Similar to the Dirac equation, the velocity operator v 
in the Kane model has a discre te sp ectrum. Each compo- 
nent of v has eigenvalues ±y/2/3V and 0, which corre- 
spond to (pure) electron, LH, and HH states. In general, 
for a wave packet containing a superposition of electron, 
LH, and HH states we have a finite probability to mea- 
sure each of these discrete values. The components of the 
velocity v do not commute. We get 



2i 



[Vj,V k \ 



■ V 2 e m 







(34a) 



where are the 4x4 spin matrices for spin s = 3/2. 
This corresponds to the uncertainty relation for j =/= k 



A Vj Av k > — 
6 



(34b) 



Note that Eq. (|34a[) implies that the minimum uncer- 
tainty depends on the dominant character of the wave 
function. The lower bound V 2 /& requires a LH state. 
For an electron state, the minimum uncertainty is 7 ?2 /3, 
whereas for a HH state it is V 2 /2. 

We omit here the lengthy expressions for r(i) and v(t). 
It follows from Eq. (|33a|) that the oscillating parts f(£) 
and v(£) of r(t) and v(t) depend on the frequencies 



E- = 2A 



K 



E±-E = E g /2 ± A K 



(35a) 
(35b) 



Unlike in the models discussed above, f 2 (t) and v 2 (t) in 
the Kane model are not diagonal in spin space. Hence 



these quantities depend explicitly on time, oscillating 
with the frequencies given in Eqs. (|35[) . However, we can 
estimate the magnitude of these quantities by neglect- 
ing the oscillatory terms and diagonalizing the resulting 
matrices. We get the following twofold-degenerate eigen- 
values for r 2 (t) 



'■(t) 




E 2 A 2 K /4 



± 



3E a 



8A K ; P 2 



A Taylor expansion shows that for small moan velocities 
( "nonrclativistic limit") we thus have two characteristic 
length scales for the oscillatory motion, the de Broglie 
wave length Xb and an effective Compton wave length 12 



X K = 



KP_ 

E n 



(37) 



We have X K ~ 7 A in GaAs and X K ~ 40 A in InSb which 
should be compared with Xd = 3.9 x 10~ 3 A. Note that, 
in the nonrelativistic limit, the de Broglie wave length 
becomes a fourfold degenerate eigenvalue of r 2 (t), i.e., 
it characterizes the oscillatory motion of electron, HH, 
and LH states. For large mean velocities ( "relativistic 
limit"), the de Broglie wave length is the only length 
scale characterizing f (£). Similarly, we get for v (t) 



v 2 (t) 



V 2 (|A 2 



K 



\E 2 g =t \EgA 



K 



(38) 



i.e., the magnitude of v is of the order of V for both 
small and large mean velocities. Again, the minimum 
uncertainty of v [Eq. (|34b[) ] and the magnitude of the 
oscillations of v are characterized by the same parameter. 
The mean velocity reads 



vi/i= \H K + ^- 



1 - 



HrE 5 



4 • (39) 



The components of v commute with each other and they 
are constants of the motion. The mean position operator 
reads 



r(i) 



vt+il 



3E„H K h 2 v 



p 2 V 2 > 4A 



K 



dHl 
dp 



2pH 2 K 



h[3A% 



(Eg/2f 



8iH 2 K A% 



(40a) 



where {A, B} = ^(AB + BA) denotes the symmetrized 



product of A and B, 
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and v denotes the acceleration v = (i/h)[HK,'v] . The 
components of f(t) do not commute with each other. We 
do not give here the lengthy expressions. 

Orbital angular momentum L = r x p and spin S also 
oscillate as a function of time. Similar to the Dirac and 
Luttinger cases, these oscillations arise even though free 
particles are considered with no external or effective mag- 
netic field present. However, the total angular momen- 
tum J = L + S does not oscillate as a function of time 
which, as always, reflects the fact that [J,Hk] = 0. 



The velocity operator, in complex notation V = v x 
is given by 



■iv v 



V(t) = - e" 1 ^* , 
m 

so that v(t) = v(i) and v 2 (t) = (uj c A c ) 2 . 
nents v x and v y do not commute, 



[Vx,Vy] = [v x (t),Vy{t)} = 



i h eB z 



(46) 

The compo- 
(47a) 



VI. LANDAU HAMILTONIAN 

There are several remarkable similarities between the 
spin-dependent Hamiltonians discussed above and the 
well-known and rather simple case of the Landau Hamil- 
tonian 2 "£. describing the cyclotron motion of 2D electrons 
in the presence of a magnetic field B z > perpendicular 
to the 2D plane. The Landau Hamiltonian is given by 



p: 



-P 2 y 



2m 



where p is the kinetic momentum with 



\Px,Py. 



= —ih eB z 



(41) 



(42) 



For the elementary charge e we use the convention e = |e| 
The time-dependent position operator 



r(t) = r 



P • / .s . P x e 2 

smlWct) -I 

muj c muj c 



[cos(w c t) - 1] (43a) 



can be written in a compact form using the complex no- 
tation R — x — iy and P — p x — ip y (Ref. l43l ). which 
highlights the analogies between the Landau Hamilto- 
nian and the models in the preceding sections. We get: 



R(t) = R 



P 1 



-iuj c t 



tUJ c 



(43b) 



where lu c = eB z /m is the cyclotron frequency. Equation 
(|43b[) shows that P/m behaves similar to the F opera- 
tors in the preceding sections^ The magnitude of the 
oscillations of R(t) is the radius A c = p/(mu) c ) of the 
cyclotron orbit. Ignoring the oscillations with frequency 
uj c , we have 



R(t) = R 



iP 

mtUc 



= C 



(44) 



independent of t, which corresponds to the center of the 
cyclotron orbit (the guiding center). The components x 
and y of R do not commute 



[x,y] = i\\ 



(45a) 



where A c = y / h/(eB z ) is the magnetic length. Equa- 
tion (I45ap can be written as an uncertainty relation 



Ax Ay > \ X 2 C 



(45b) 



which corresponds to the uncertainty relation 
Av x Av y > | (lo c A c ) 2 , 



(47b) 



which should be compared with Eqs. fTS]) and (f2"3"]) . Obvi- 
ously, implications arising from this uncertainty relation 
become relevant only for sufficiently large magnetic fields 
when A c becomes comparable to A c . 

The velocity V(t) is not a conserved quantity, which 
reflects the effect of the Lorentz force. We have 



V(t) 



-iui r P 



The mean velocity operator vanishes, 
V = , 



(48) 



(49) 



because, on average, the particle is at rest for B z ^ 0. 
This also implies [v x ,v y ] = 0. 

Our analysis indicates that the dynamical properties 
of the Landau model bear strong resemblances to those 
exhibited by models showing zitterbewegung-like motion. 



VII. LANDAU-RASHBA HAMILTONIAN 

An interesting example combining two types of oscil- 
latory motion can be found by considering the interplay 
between 2D cyclotron motion (Sec. IVI|) and Rashba spin 
splitting (Sec. IIII|) . The Hamiltonian for that situation 
reads 



H, 



cR 



H r + Hi 



Mb 



cr?B, 



(50) 



Here we have also included a Zeeman term with Lande 
factor g and Bohr magnetic moment /ib = eh/(2m e ) 
(where m c denotes the electron mass in vacuum), and 
the terms Br and H c are given in Eqs. (|19bp and (|41jl . 
For the following calculation we replace the components 
p x and p y of the kinetic momentum by creation and an- 
nihilation operation operators for Landau levels, a) and 
a, defined in the usual way, 



Ac£ 

\/2h 



(51) 



and a) is the adjoint of a. The resulting expression 
for B c n (Ref. [27l ) is equivalent to the Jaynes-Cummings 
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model 4 ^ in the rotating-wave approximation. To find the 
time evolution of the observables in the Heisenberg pic- 
ture, we first separate H c r into two commuting parts, 



l cR 



H^2, where 



H 



(i) 

cR 



= huj c ( a) a + 



H 



(2) _ 
cR ~ 



V2ih( 



1 + ct, 



-aV_)-^(, 



gm 
2m r 



(52a) 



(52b) 



Here we used a± = (a x ± ia y )/2. 

It is straightforward to calculate the time evolution of 
the spin component parallel to the magnetic field, 



S z (t) = S z - icr z H R 



1 



(53) 



This result is the generalization of Eq. (I24b[) to the case 
of a finite magnetic field. Interestingly, time averaging 
the r.h.s of Eq. (|53j) does not result in a vanishing spin 
component parallel to the field direction. We find 



h 2 u c 



AH 



(2) 
cR 



1 



gm 
2m r 



(54a) 



Neglecting Zeeman splitting and considering the limit of 
small B-,, this result becomes 



S z 



Aap 2 



(er x p) • e 2 



(54b) 



which is exactly the finite value of the spin component 
parallel to the magnetic field that was obtained in semi- 
classical calculations of spin-split cyclotron orbits4£ 

To calculate the time evolution of the position opera- 
tor, we use the complex notation from Sec. IVIl We have 



R = C 



iP 



muj r 



= C + iV2\ c a 



(55) 



where C is the position of the guiding center, see Eq. 
(|44|) . Even in the presence of Hr, the guiding center C 
remains a constant of the motion, [C, H c r] — 0. The 

time evolution of P due to is just a trivial factor 
e _I " ct , so that we only need to evaluate the time evolu- 
tion ofFoca under Hf R \ (Note that [H^ R >, H%>] = 0.) 
This problem has been solved for the Jaynes-Cummings 
model . 45 ' 47 Translating into our situation, we get for the 
time-dependent position operator 



R(t) = C 



i exp | 



-i( We + W+ )t] f^P +2iaa _ 



i exp [— i(uj c + t-<J-)i\ 



h 2ia<7- 

ijj c m 



(56a) 



with the frequency operators lu± given by 



hu ± = -H™ ± V {H%)" + 2Tuo c ma 



(2) A 2 



(56b) 



The terms proportional to er_ in Eq. (|56a|) are reminis- 
cent of the oscillatory motion in the Rashba case for 
B z = 0, where the amplitude of the oscillations is in- 
versely proportional to the de Broglie wave length and 
independent of a, see Eq. (|22[) . Here these terms con- 
tribute to a spin-dependent renormalization of the cy- 
clotron radius. We also note that R(t) = C, so that Eq. 
(|45| remains valid in the presence of Hr. 
The velocity operator is given by 



P 

V = R = 2iaa- . 

m 

The commutator of the components of V, 
. . iheB 



(57) 



[Vx,Vyl 



11) ~ 



+ 2ia 2 a z , (58a) 



is the sum of the corresponding results obtained sepa- 
rately from H c and Hr [see Eqs. (|23b[) and (|47a[) ]. How- 



ever, in the uncertainty relation 

Av x Av y > 

" 2m 2 



a 



(58b) 



the two contributions are subtracted, thus reducing the 
minimum uncertainty. The time dependence of V can 
be readily obtained by taking the time derivative of Eq. 
(|56ap . It can be written as 



V(t) = e _i ^ c+w+ ' t F + + e ~ i ( u '<' +u) -) t F_ 
with the complex amplitude operators 



(59a) 



V ( H [2 R ] +2ma 2 J'' 



LU+ — U>- 



2ia H %-^ 

LU+ — U>- 



(59b) 



VIII. CONCLUSIONS AND OUTLOOK 

We studied a variety of qualitatively different model 
Hamiltonians for quasi-free electrons that exhibit zitter- 
bewegung-like oscillatory motion. A number of features 
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can be identified that are widely shared as discussed in 
Sec. HI Here we finally point out open questions. 

For the Dirac Hamiltonian, the amplitude of the zit- 
terbewegung of r(t) is given by the Compton wavelength 
in the nonrelativistic limit and by the de Broglie wave 
length in the relativistic limit. For those Hamiltonians 
having a gap that vanishes for p — > 0, the length scale 
of oscillations in r(t) is always given by the de Broglie 
wavelength Xb = h/p, independent of the magnitude of 
spin-orbit coupling. It is surprising that the amplitude 
of the oscillations of r(t) diverges in the nonrelativistic 
limit p — > 0. 

The most interesting but also, at least in our present 
work, a largely open aspect is the experimental observ- 
ability of zitterbewegung-like effects. Certainly, any mea- 
surement of the oscillatory motion must obey the fun- 
damental uncertainty relations [Eq. Q] discussed in our 
work. Furthermore, we have already commented on the 
intimate relation between oscillations in position and spin 
space. However, while spin precession due to spin-orbit 
coupling can be observed experimentally, 25 it is often ar- 
gued that the zitterbewegung of r(i) "is not an observ- 
able motion, for any attempt to determine the position of 
the electron to better than a Compton wavelength must 
defeat its purpose by the creation of electron-positron 
pairs" (Ref . @) . We note that the same argument can be 
applied to Bloch electrons in solids where electron-hole 
pairs can be created*^ 
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APPENDIX: IMPORTANT FORMULAE 

Here we briefly summarize important formulae that 
are used in our discussion of the oscillatory motion in 
various models. The Heisenberg equation of motion for 
an operator A reads 

It has the formal solution 

A(t) = e lHt / n A(0) e' lHt/h . (A.2) 

In particular, the velocity operator v is defined by the 
Heisenberg equation of motion for the position opera- 
tor r, 

Throughout we use the convention that A(t) denotes an 
operator in the Heisenberg picture and A = A(0) is the 
corresponding operator in the Schrodinger picture. 

In general, the uncertainty principle for two noncom- 
muting observables A and B reads 

AAAB > \ |([A,B])| , (A.4) 

where the uncertainty AA of A is defined as 

AA = v/(A 2 ) - (A) 2 . (A.5) 
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